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ABSTRACT 

We combine DensePak integral field unit and TAURUS Fabry-Perot observations of 13 nuclear 
rings to show an interconnection between the kinematic properties of the rings and their resonant 
origin. The nuclear rings have regular and symmetric kinematics, and lack strong non-circular 
motions. This symmetry, coupled with a direct relationship between the position angles and 
ellipticities of the rings and those of their host galaxies, indicate the rings are in the same plane 
as the disc and are circular. From the rotation curves derived, we have estimated the compactness 
(u 1 2 /r) up to the turnover radius, which is where the nuclear rings reside. We find that there is 
evidence of a correlation between compactness and ring width and size. Radially wide rings are 
less compact, and thus have lower mass concentration. The compactness increases as the ring 
width decreases. We also find that the nuclear ring size is dependent on the bar strength, with 
weaker bars allowing rings of any size to form. 

Subject headings: Galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure - galaxies: 
nuclei 


1. Introduction 

Nuclear rings in the central kiloparsecs of galax- 
ies are usually observed as conglomerations of 
young massive stars in distinct compact groupings 
(Buta & Combes 1996; Knapen 2005; Comeron 
et al. 2010). The star forming rings are believed 
to form as a result of radial gas inflow towards 
the central region, which stagnates near dynami- 
cal resonances. In barred galaxies the inflow orig- 
inates in the gravitational torques set up by the 
bar (e.g., Athanassoula (1992); Heller & Shlos- 
man (1994, 1996); Piner et al. (1995); Knapen et 


al. (1995); Buta & Combes (1996); Regan et al. 
(1997); Benedict et al. (2002)). Some 20% of nu- 
clear rings occur in non-barred galaxies (Knapen 
2005; Comeron et al. 2010) and in those cases the 
inflow and resonances may well be related to the 
presence of weak ovals, past interactions or merg- 
ers, or even strong spiral arms. A resonant origin 
of all nuclear rings is most plausible. 

The radial location of a nuclear ring can be de- 
scribed by the position of the inner Lindblad res- 
onance^) (ILRs) (Sirnkin et al. 1980; Shlosman 
et al. 1989; Buta & Crocker 1993; Friedli & Benz 
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1993; Heller & Shlosman 1994; Knapen et al. 1995; 
Buta & Combes 1996), which is typically at ~1 
kpc from the galaxy center although smaller and 
larger nuclear rings exist (Buta & Crocker 1993; 
Comeron et al. 2010). In a non-axisynnnetric po- 
tential, the nuclear ring forms near the ILR(s), 
where the inflowing molecular gas slows down, 
accumulates, and initiates star formation. For 
barred spiral galaxies, the ring can be formed in- 
side the large scale galactic spiral structure if the 
bar pattern speed, fib, closely matches the inner 
Lindblad precession frequency (ft — ac/2), where 
fl is angular velocity and k is the radial epicyclic 
frequency (Athanassoula 1992). 

Only a few in-depth observational studies of the 
kinematics of nuclear rings at high angular resolu- 
tion exist, but these have revealed some consistent 
kinematic characteristics of nuclear rings. Bene- 
dict et al. (1996) showed that non-circular motions 
near the intersection of the ring and bar dust lanes 
correspond to radially inflowing molecular gas that 
is apparently feeding the ring in NGC 4314. Rey- 
naud & Downes (1998) concluded for NGC 1530 
the existence of weak velocity transitions at the 
ends of the bar which are inversely correlated with 
strong star formation in the same areas, including 
the nuclear ring location. Regan et al. (1997) fur- 
ther note an increase in the residual velocity along 
an axis perpendicular to the bar. Knapen et al. 
(2000) (see also Knapen et al. (1995)) and Zu- 
rita et al. (2004) used CO and Ha emission in the 
central region of NGC 4321 and NGC 1530, re- 
spectively, to reveal gas streaming inward towards 
the ring along the bar dust lanes. Jogee et al. 
(2002) studied the fueling mechanisms of the ring 
in NGC 5248 to deduce that large star forming 
clusters in the ring have been triggered by a bar- 
driven spiral density wave. More recently, Allard 
et al. (2006) follow up on the ring in NGC 4321 to 
show relatively low H j3 gas dispersion within the 
ring, which can be an indicator of active star for- 
mation and cool gas inflow from which the massive 
stars inherited the low velocity dispersion. 

The existing kinematic studies provide insight 
into a small number of nuclear rings, but fail to 
reach more general conclusions applicable to the 
overall population. Apart from the differing char- 
acteristics of the data used in these studies, this is 
due mostly to the intrinsically small angular size 
of nuclear rings, and the dusty environments in 


which they exist, both of which make rings non- 
trivial to observe with most common instruments 
yielding kinematic data. In this paper we endeavor 
to start rectifying this situation by presenting a 
more homogeneous data set on thirteen nuclear 
rings. 

The nuclear rings studied spectroscopically in 
the current paper were all part of the sample of 
Mazzuca et al. (2008) (hereafter known as M08), 
who provided photometric insight into the strong 
star forming nature of nuclear rings. From their 
Ha observational study of 22 nuclear rings, they 
find that many of the rings exhibit a well-defined 
age distribution pattern. M08 conclude that this 
occurrence can be a result of the combination of 
bar-induced dynamics and gravitational instabil- 
ities which are occurring in the proximity of the 
rings. These results add credence to the idea that 
the mass inflow driven by the bar along the dust 
lanes to the ring contact points is a key require- 
ment for the fueling of active (i.e., star forming) 
nuclear rings (and might also yield a possible sig- 
nature in the form of non-circular motions). Fur- 
ther evidence for this idea from detailed observa- 
tions of a handful of galaxies has been presented 
by, e.g., Allard et al. (2006); Boker et al. (2008), 
and Knapen et al. (2010). 

We now use kinematic information within and 
just outside the nuclear rings in 13 galaxies to in- 
vestigate the possible relations between the optical 
morphology of the rings, the parameters derived 
from the rotation curve, and their resonant na- 
ture. We begin by presenting an overview of the 
observations and morphological properties of the 
Ha sample in Section 2. We discuss the basic data 
reduction method, which includes the construc- 
tion of the velocity field maps, in Section 3. In 
Section 4 we discuss the kinematic parameters as- 
sociated with the ring rotation curves, and present 
results pertaining to non-circular motions in and 
around the nuclear rings, including velocity resid- 
uals. Our detailed analysis of the results is dis- 
cussed in Section 5, with concluding remarks in 
Section 6. 

2. Sample and Observations 

The galaxies studied here (see Table 1) form 
a subset of the sample presented in the imaging 
survey of M08 (which was based on B, /, and 
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Ha imaging). The data in M08 were obtained to 
characterize the morphology of the nuclear rings, 
and to analyze the age distribution of the indi- 
vidual stellar clusters forming each ring using the 
Ha equivalent width. Results from this photo- 
metrically based survey reveal that the rings con- 
tain very young (1 Myr to 10 Myrs) hot massive 
stars, with some of the youngest clusters located in 
proximity to one or both of the intersection points 
of the bar dust lanes and the outer edge of the 
ring. The ring ellipticity and position angle ap- 
proximately match that of the respective galac- 
tic disk, indicating that the rings are in the same 
plane as the disk and are circular (M08; see also 
Comeron et al. (2010)). 

The kinematic sample was chosen from the 
imaging parent sample based on unambiguous de- 
tection of H II regions delineating the ring, a non- 
zero inclination (z), and a measurable ring posi- 
tion angle (<£ r ). Rings too diffuse or small, or 
with no resolved Hu regions detected using the 
method described in M08 (i.e. , weak signal and 
resolution), were rejected. With the Densepak in- 
strument we cannot detect any Hil regions for 
nuclear rings with angular sizes less than 0"5 in 
radius, however, we have previously unpublished 
high resolution Fabry-Perot Ha data of NGC 1300 
and NGC 6951 from the TAURUS instrument on 
the 4.2 m William Herschel Telescope, which we 
include in this sample. Similarly, we have used 
the TAURUS data for NGC 4321 from Knapen 
et al. (2000), a nuclear ring which meets our cri- 
teria but was added to the M08 sample after the 
DensePak observing runs. Lastly, although not in- 
cluded in the M08 sample, we observed two well 
known nuclear rings (NGC 2903 and NGC 3351) 
during the December 2003 run, which we also add 
to this sample. 

We obtained two-dimensional Ha velocity fields 
of the sample over two observing runs, in De- 
cember 2003 and April 2004, using the DensePak 
fiber-optic array and bench spectrograph on the 
3.5 m WIYN 1 telescope. The DensePak integral 
field unit (IFU) spectrograph (Barden & Wade 
1988) consists of 91 fibers bonded into a 7 x 13 
staggered rectangular grid that covers an area of 
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30 x 45arcsec of sky with center-to-center fiber 
(each fiber is 2.82arcsec in diameter) spacings of 
3.75arcsec. Four of the fibers are allocated as sky 
fibers, and are spaced around the grid roughly 1' 
from the center. The spectrograph was configured 
with the Bench camera, which uses a Tek 2048 
CCD (T2KC), with the 860 (mm -1 grating, in 
the second order, centered at 6575 A. We used a 
third order blocking filter (GG495). This arr ange- 
ment provided a final spectral coverage of 925 A, 
from 6085 A to 7010 A, at a resolution of 0.45 A per 
pixel. The spectral range included Ha as well as 
the [N n] AA6548, 6583 and [S n]AA6716, 6731 lines. 
The effective instrumental resolution was different 
for the two runs, which is likely due to differences 
in the instrument setup. For the December 2003 
run, the instrumental resolution (FWHM) aver- 
aged 0.97 A (44 kms -1 ) around Ha, with the res- 
olution for the April 2004 run 1.31 A (60 kms -1 ). 

The array contains six dead fibers which de- 
creases the number of active fibers to 85. To avoid 
these areas, the projected image of each nuclear 
ring was placed in the upper half of the fiber ar- 
ray. Because of this method, many of the rings 
were close to the edge of the field of view (FOV) 
along the array minor axis, which limits our dis- 
cussion of the environment exterior to the rings. 
We used two to four pointings, each of 1800 s in 
duration, offset to fill the area of the dead fibers. 
This method also improved the sampling in some 
cases. For all exposures, the major axis of the 
array was aligned along a North-South direction. 

The Fabry-Perot observations for NGC 4321 are 
detailed in Knapen et al. (2000), with very sim- 
ilar procedures implemented for NGC 1300 and 
NGC 6951. The pixel size in the reduced Fabry- 
Perot cubes is 0.27 arcsec squared, and the sep- 
aration between velocity channels is 15.7 kms -1 . 
Each cube has 55 planes, thus covering a total ve- 
locity range of just over 1000 kms -1 . Continuum 
subtraction was performed using the channels on 
either extreme end of this range. 

3. Data Reduction 

The raw CCD frames of the DensePak data 
were zero- and bias-corrected, flat-fielded, cosmic- 
ray cleaned, and combined using the standard 
iraf no ao packages (zerocombine, ccdproc, 
FLATCOMBINE, COSMICRAYS, and IMCOMBINE, re- 
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Table 1: Morphological and photometric parameters of DensePak sample 


NGC 

Morph 

Offset 

i 

4>r 

<t>d 

4>b 

e r 


Radius 

Radius 

Run 


Type 

X 

y 









Date 



n 

(") 

(°) 

(°) 

(°) 

(°) 



(") 

(kpc) 


(1) 

(2) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

473 

SAB(r)0/a 

-4.4 

+7.2 

49 

154 

153 

164 

0.37 

0.37 

12.2 x 6.9 

1.7 x 1.0 

Dec03 

1300 

SB(rs)b 

n/a 

n/a 

40 

135 

106 

102 

0.25 

0.34 

4.1 x 3.1 

0.3 x 0.2 

Sep98 

1343 

SAB(s)b 

-4.4 

+7.8 

41 

60 

80 

82 

0.25 

0.37 

8.8 x 6.6 

1.2 x 0.9 

Dec03 

1530 

SB(rs)b 

-1.6 

+11.3 

48 

25 

8 

122 

0.35 

0.29 

6.8 x 4.9 

1.2 x 0.8 

Dec03 

2903 

SAB(rs)bc 

-6.8 

+8.3 

65 

8 

17 

24 

0.4 

0.52 

1.7 x 1.0 

0.06 x 0.04 

Dec03 

3351 

SB(r)b 

-4.2 

+6.1 

34 

20 

13 

112 

0.3 

0.32 

2.7 x 1.9 

0.15 x 0.11 

Dec03 

4303 

SAB(rs)bc 

-3.3 

+7.5 

18 

88 

- 

10 

0.14 

0.11 

3.3 x 2.8 

0.2 x 0.2 

Apr04 

4314 

SB(rs)a 

-3.0 

+6.3 

26 

135 

- 

135 

0.1 

0.11 

6.6 x 5.9 

0.3 x 0.3 

Apr04 

4321 

SAB(s)bc 

n/a 

n/a 

30 

170 

30 

153 

0.12 

0.15 

8.8 x 7.0 

0.7 x 0.6 

May95 

5248 

SAB(rs)bc 

-3.3 

+6.5 

46 

115 

110 

137 

0.3 

0.28 

6.6 x 4.6 

0.7 x 0.5 

Apr04 

5953 

SAa 

-6.0 

+6.6 

26 

172 

169 

no bar 

0.1 

0.17 

6.1 x 5.5 

1.0 x 0.9 

Apr04 

6951 

SAB(rs)bc 

n/a 

n/a 

45 

146 

170 

85 

0.2 

0.17 

4.6 x 3.7 

0.5 x 0.4 

Sep93 

7742 

SA(r)b 

+0.7 

+7.5 

18 

133 

- 

no bar 

0.05 

0 

9.9 x 9.4 

1.0 x 1.0 

Dec03 


Notes: Morphological and photometric parameters for the observed sample. Galaxies are listed by NGC 
number in order of increasing RA (col. 1) with morphological type (col. 2) from de Vaucouleurs et al. (1991; 
hereafter RC3). Col. 3 lists the offset of the optical center of the Ha image from the center of the DensePak 
array; this is not applicable to the TAURUS data for NGC 1300, NGC 4321, and NGC 6951. The nuclear 
ring inclination i (col. 4) and photometric ring position angle fa (col. 5) have been derived from Ha imaging 
by M08). The disk position angle fa (col. 6) is from the RC3 with a dash indicating that a ring is circular 
with no definable position angle. The bar position angle fa (col. 7) and photometric ring ellipticity e r (col. 8) 
are from M08. Disk ellipticity ea (col. 9) is from RC3. The radius of the ring, shown as the semi-major axis 
by the semi-minor axis (col. 10 and col. 11), was derived from the imaging data of M08 using the distance 
given there. The ring semi-major axis and ring ellipticity values are all consistent with those reported in 
Comeron et al. (2010), taking into account that the methodology applied was slightly different. Run dates 
are in col. 12. 


4 



spectively). The combined images were wavelength- 
calibrated using the NOAO iraf data reduction 
package dohydra to extract the spectrum yielded 
by each fiber. A CuAr comparison source was 
used for wavelength calibration. The reduced 
two-dimensional spectral image contains 89 rows 
(comprising 85 data fibers and 4 sky fibers), where 
each row is one pixel in height and corresponds to 
a single fiber. 

We created GIPSY 2 scripts to identify and then 
subtract the continuum and sky from the spec- 
trum following the method described in Bershady 
et al. (2005). Due to pointing uncertainties in- 
herent to the instrument, we established the mea- 
sured pointing offsets (i.e., the shift from the cen- 
ter of the DensePak fiber array to the ring center) 
from the observations themselves. We followed 
the procedure outlined in Swaters et al. (2003), 
in which the continuum levels measured from the 
spectra are compared with levels derived from a 
reference image with known coordinates. For our 
sample we used the /-band images from M08. By 
repeating this process for all positions near the 
nominal pointing, we created a grid of y 2 values, 
in which the minimum y 2 corresponds to the mea- 
sured pointing offset. We estimate that the un- 
certainty on the derived pointing is approximately 
0.5 arcsec. Because the sky fibers contained no Ha 
emission from the target galaxy, we used the aver- 
age of the four reserved DensePak sky fibers, each 
located 1 arcmin from the center of the array. 

To combine the pointings into a sparse velocity 
field, we first fitted Gaussians to the continuum- 
and sky-subtracted Ha emission line profiles, using 
a 3cr threshold to avoid spurious detections. The 
fitted velocities derived for each fiber were placed 
in a map at corresponding fiber positions. The 
observed radial velocities were corrected to helio- 
centric velocities. In the cases of multiple velocity 
measurements for the same pixel in the velocity 
field, we averaged the measured velocities. To aid 
in the visual representation of the sparse velocity 
field, we also constructed interpolated smoothed 
velocity field maps. We replaced each data point 
in the sparse velocity field with a Gaussian beam 
of 5.5 arcsec and weighted the overlapping sections 
between the points by the relative intensities of the 


2 Groningen Image Processing System, 
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overlapping Gaussians. The resulting smoothed 
velocity fields are shown in Figure 1 (left, panel 1). 
Because of interpolation effects these contiguous 
velocity fields may be uncertain, especially near 
the edges. Any quantitative interpretation comes 
from the sparse velocity maps, which are not in- 
terpolated. 

The Fabry-Perot data were wavelength and 
phase calibrated, with sky and continuum emis- 
sion subtraction performed, as described in detail 
in Knapen et al. (2000). The resulting reduced 
data cubes were used to produce Ha velocity field 
maps using GIPSY, and following the prescriptions 
given by Knapen (1997) and Knapen et al. (2000). 
The resulting velocity fields are shown in Figure 2. 

4. Analysis 

4.1. Rotation Curves 

Rotation curves are powerful tools for revealing 
intrinsic kinematic characteristics of galaxies, in- 
cluding those of the nuclear rings they may host. 
In our sample, the range in radius spanned by 
the rotation curves is small when compared to the 
overall disk size, but this is a result of the small 
field of view of the instrument and of the sensitiv- 
ity limit. 

For each galaxy we derived the rotation curve 
(Figure 3) by fitting tilted rings to the sparse ve- 
locity field through an iterative process (see Bege- 
man (1989) for procedural details). We chose 
the incremental radius for each concentric ring to 
be 1.5 arcsec, which is approximately half of the 
width of a given fiber. We supplied the initial in- 
put values for the nuclear ring center, inclination, 
and position angle originate from /-band images 
using the AuxPort camera as observed in M08. 
We derived the initial center position by taking 
the point where the /-band emission is the highest 
near the galaxy center. The input value for v sys in 
the fit to the velocity field is from the NASA/IPAC 
extragalactic database (NED). One exception is 
NGC 5953, whose systemic velocity was derived 
from the data due to a 26 kms -1 difference with 
the NED value. For the initial run, although we 
supplied input values, all parameters were chosen 
to be free. We first determined the center posi- 
tion and then fixed only that value for the second 
iteration to then determine the systemic velocity. 
Once the center and systemic velocity were fixed 
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Fig. 1. — Results from the DensePak sample. From left to right, I band image from M08 with velocity field 
overlaid in contours, residual velocity field with contours, and Ha emission image from M08. The center of 
the galaxy is indicated with a cross in each panel. Contour levels in the left panel are spaced by 5 kms -1 
with dotted contours representing velocities less than the systemic velocity and solid contours representing 
those velocities greater than the systemic velocity. Residual velocity contours and grayscales within the rings 
generally range from —10 to 10 kms -1 except for NGC 1530 whose contours range from —20 to 20 kms -1 , 
and NGC 473 and NGC 7742 which both have residual ranges from —5 to 5 kms -1 . N is up, E to the left. 
The size of the field shown is lOarcesc (RA) by 40arcsec (DEC) for all galaxies. 


6 














Fig. 2. — As Fig.l, now for the three galaxies observed with the TAURUS Fabry-Perot instrument. Panels 
are, from left to right, velocity field (grayscale and contours), result after subtracting a model velocity field 
based on the rotation curve from the observed velocity field, and Ha integrated intensity. The center of 
the galaxy is indicated with a cross. The field of view shown is 15, 35, and 15arcsec square for NGC 1300, 
NGC 4321, and NGC 6951, respectively. Contour levels for the first (leftmost) panel are for NGC 1300 from 
1472 to 1552kms _1 (white contours) and from 1572 to 1672kms _1 (black); for NGC 4321 from 1480 to 
1560 (white) and from 1580 to 1680kms _1 ; and for NGC 6951 from 1300 to 1440 (white) and from 1460 
to 1600kms _1 —all contours are spaced by 20kms _1 and the grayscale range is the same as that spanned 
by the contours. In panel 2, contour levels for all galaxies are from —45 to Okms -1 (white) and from 15 to 
45kms _1 in steps of 15kms _1 ; the grayscale range is from —50 to 50kms -1 for NGC 1300 and NGC 6951 
but from —30 to 30kms _1 for NGC 4321. For all three galaxies, the grayscale in panel 3 indicates the range 
of values in instrumental units. 
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we determined the inclination and position angle 
(third and fourth iterations, respectively). With 
the four values then fixed, we fit the rotational 
velocity. We gave more weight to positions along 
the major axis, since those data points carry most 
of the information about rotational velocities. In 
most cases we used the inclination of the outer disk 
for the final fit. However, in the cases of low incli- 
nation (below 25°), we adopted the photometric 
orientation parameters from M08 as initial inputs 
to the routine since the fitting routine produced 
inclinations with high uncertainties. Even with 
such low inclinations the position angle was deter- 
mined accurately with a tilted ring fit. Because 
our data sets span a small range in radius, it was 
not feasible to make sophisticated fits to the ve- 
locity fields. 

The uncertainties on the velocities derived from 
the line fits range from 0.5 kms” 1 to 5 kins -1 . An 
additional source of uncertainty is due to the un- 
known location of the Halpha emission within the 
area seen by each fiber. The amplitude of this un- 
certainty depends on the velocity gradient across 
the face of the fiber. We estimated this gradient 
for all fibers using the derived rotation curves, and 
found that a value of 8 kms” 1 is representative of 
a ler uncertainty. 

Rotation curves for the Fabry-Perot data were 
similarly generated, but here, the smaller pixel size 
and full spatial sampling allowed us to fit the ro- 
tation curves with radial steps of 0.5arcsec from 
the center outward. 

4.2. Residual Velocity Fields 

We constructed residual velocity maps (Fig- 
ure 1, panel 2) for the velocity fields (both the 
sparse and interpolated versions for the DensePak 
data) by subtracting a smooth model velocity 
field, constructed from the rotation curve and the 
assumed orientation parameters. As with the in- 
terpolated velocity fields, the interpolated residual 
maps should only be used for general impressions. 
These maps, however, are useful because they give 
a better sense of patterns in the image. 

4-2.1. Within the rings 

We find that all of the nuclear rings in our sam- 
ple have low residual velocities (both approaching 
and receding) near the central ridge of the ring 


(ranging from ~5 km s 1 to ~20kms - 1 ). The nu- 
clear rings in NGC 473 and NGC 7742 are not 
resolved enough to comment on here, as the asso- 
ciated residuals do not overcome the uncertainties 
associated with the measured velocity, wavelength 
calibration, and position of the Ha peak within a 
given fiber beam. Patches of approaching and re- 
ceding residual velocity peaks of ~10kins _1 occur 
in the nuclear rings of NGC 1343, NGC 4314, and 
NGC 5953, and reach to near 20 kms” 1 in the 
rings of NGC 1530 and NGC 5248. Since these 
values are close to the combined measurement un- 
certainty, they can provide no firm evidence for 
the existence and origins of non-circular motions 
(in NGC 4321 the situation is better, as discussed 
in detail by Knapen et al. 2000). 

A comparison of the kinematic and photometric 
centers of the sample reveals no strong asymme- 
tries. The degree to which the centers differ is pro- 
portional to the lopsidedness of the potential, and 
thus a good indicator of the amount of non-circular 
motion (Franx et al. 1994). After computing the 
difference between the kinematic and photometric 
center positions (see Table 2), we see that the cen- 
ter positions agree to within the uncertainties of 
the photometric center measurements. This agree- 
ment, in combination with the small noncircular 
motions, suggests the nuclear rings are circular in 
nature. 

4-2.2. In the vicinity of the rings 

For some of the galaxies in our sample, our ob- 
servations cover part of the velocity field outside of 
the nuclear ring as well. This allows us to study 
the interaction between the bar and the nuclear 
ring in those cases. 

Sufficient, albeit limited, Ha radial coverage al- 
lows us to comment on the velocities of the barred 
galaxies near the outer edge of the nuclear ring, 
perpendicular and parallel to the bar major axis. 
In all cases, we can see the influence of the strong 
bars as they create high residual velocities. How- 
ever, once they interact with the nuclear ring outer 
boundary, residual velocities sharply decrease and 
plateau within the nuclear rings. For NGC 4314, 
we see ordered velocity patterns increasing near 
the outer edge of the ring (both approaching and 
receding) close to both sides of the bar minor axis, 
where the bar dust lanes merge with the ring. 
Residuals approach 60 kms” 1 in both directions, 
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Rotational Velocity (km/s) Rotational Velocity (km/s) Rotational Velocity (km/s) Rotational Velocity (km/s) 


N473 


N 1 300 


N 1 343 




N 1 530 


N2903 


N3351 




N4303 


Radius (kpc) 
N4314 


Radius (kpc) 

N4321 





0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 

Radius (kpc) Radius (kpc) Radius (kpc) Radius (kpc) 


Fig. 3. — Rotation curves with the radial range (i.e. , width) of the nuclear rings shadowed. The turnover 
radius is denoted by the dashed vertical line. See Table 3 for ring widths. 
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Table 2: Kinematic parameters for the observed sample 


NGC 

Kiys 

Offset 

/'kin 

| 0kin ^phot | 

| <t>d ~ </kin 

| Ctl' kin 

Ctlphot | 



X 

y 




X 

y 


(kms -1 ) 

(") 

(") 

(°) 

(°) 

(°) 

(") 

n 

(1) 

(2) 


(3) 

(4) 

(5) 

(6) 


(7) 

473 

2135 

-5.0 

+7.5 

156 

2 

3 

0.6 

0.3 

1300 

505 

n/a 

n/a 

80 

55 

26 

n/a 

n/a 

1343 

2215 

-5.3 

+6.5 

7 

53 

73 

0.9 

1.3 

1530 

2450 

-1.6 

+11.0 

5 

20 

3 

0 

0.3 

2903 

560 

-4.0 

+7.4 

5 

3 

12 

2.8 

0.9 

3351 

764 

-4.2 

+7.0 

30 

10 

17 

0 

0.9 

4303 

1560 

-3.4 

+ 7.5 

117 

29 

— 

0.1 

0 

4314 

979 

-4.0 

+5.5 

122 

13 

— 

1.0 

0.8 

4321 

1578 

n/a 

n/a 

150 

20 

120 

n/a 

n/a 

5248 

1149 

-2.9 

+6.3 

115 

0 

5 

0.4 

0.2 

5953 

1991 

-6.0 

+6.7 

49 

57 

120 

0 

0.1 

6951 

1423 

n/a 

n/a 

136 

10 

106 

n/a 

n/a 

7742 

1658 

+0.7 

+6.6 

131 

2 

- 

0 

0.9 


Notes: NGC numbers of the target galaxies are listed in order of increasing RA (col. 1). Systemic velocity 
(col. 2), nuclear ring kinematic center position with respect to the center of the DensePak array (col. 3) 
and kinematic position angle (col. 4) were derived using the GIPSY ROTCUR task. Typical uncertainties for 
14 ys are ±8kms _1 , center position ±0.09arcsec, and position angle ±5°. The difference (absolute value) 
between the kinematic and photometric position angles (see Table 1) is in col. 5, with the difference between 
kinematic and disk position angles in col. 6. The difference between kinematic and photometric centers is in 
col. 7. 
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with a strong velocity transition zone near the 
outer boundary of the nuclear ring, where values 
reach 90 kms -1 (approaching) at the edge of the 
FOV. Although the spatial sampling is poor at the 
image edge, our observations do agree with veloc- 
ity values from Benedict et al. (1996). 

The residuals for NGC 5248 near the western 
side of the bar minor axis reach 40 kms -1 , with a 
weaker-defined, although evident, transition zone 
as radii approach the outer ring boundary. These 
results are consistent with the two-dimensional 
CO velocity field presented by Jogee et al. (2002). 
In both galaxies, good agreement with the litera- 
ture adds confidence that the results we see at the 
field edges are not artificial, even though our data 
are spatially sparse. 

The residual map of NGC 1530 provides us with 
the clearest picture of the interaction between the 
bar and the ring. In Fig. 4 we see large velocity 
excesses (~110kms -1 ) to the west of the nuclear 
ring, near its exterior edge. The high residuals are 
associated with the large scale bar. The strong 
residuals sharply decrease to 50km s -1 over a small 
transition zone near the outer boundary of the 
nuclear ring. Although we see no obvious resid- 
ual patterns near the bar major or minor axes of 
the nuclear ring, we note that the largest residual 
velocities (~120kms -1 ) across our sample occur 
in this galaxy, between the bar major and minor 
axes, on the west side just exterior to the ring 
outer boundary. A closer look in the vicinity radi- 
ally shows a very sharp transition at the location 
of the nuclear ring, where residual velocities drop 
to 10 kms -1 . Zurita et al. (2004) observe a nearly 
identical pattern in the corresponding Ha map of 
the circumnuclear environment. 

The TAURUS data of NGC 4321 used here were 
analyzed and discussed in detail by Knapen et al. 
(2000). In the residual velocity map, they high- 
lighted the contributions from the inner part of 
the bar and the spiral arm fragment coming into 
the nuclear (pseudo-) ring from the main part of 
the bar. These result in a broadly symmetric pat- 
tern of positive and negative residual velocities, 
mapping out the non-circular motions due to bar 
and spiral streaming motions. 



within (top plot) and exterior to (bottom plot) 
the nuclear ring. The radial coverage within the 
nuclear ring is further divided to show velocities 
along the inner and outer halves of the ring. Out- 
side of the nuclear ring, velocities from the nu- 
cleus to the ring inner edge are plotted, as well 
as those from the ring outer edge to the radial 
limits of the FOV. The velocity is plotted with re- 
spect to the systemic velocity, normalized by the 
rotation velocity, and adjusted by removing the 
inclination in order to allow comparison with the 
curve representing circular motion. The bar ma- 
jor axis (solid lines) and minor axis (dotted lines) 
are drawn for context. Note that there is a phase 
shift of 90°, plus the difference between the kine- 
matic position angle of the velocity data points 
and the photometric position angle used to cre- 
ate the model curve, to allow the cross-over point 
to occur at 180°. The area azinruthally located b 
etween the bar major and minor axes shows the 
largest velocities of the sample, with deviations 
reaching 110 km s - 1 . 
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5. Results 

5.1. Circular Nuclear Rings 

For most of the sample, we find that the nuclear 
rings are nearly circular and in the same plane as 
the disk. The shape and orientation with respect 
to the host disk were deduced by comparing the 
position angle (of the nuclear ring, <p r and disk, 
(f>d) combined with the ellipticity of each nuclear 
ring, e r to that of its host disk, We used the 
photometric position angles and ellipticities of the 
nuclear rings from M08. The kinematic position 
angles of the nuclear rings were extracted from 
the rotation curve generation. In Figure 5, we 
compare the various position angles and see that 
in most cases the photometric and kinematic po- 
sition angles are similar. The same can be said 
for both the photometric and kinematic position 
angles with respect to the galactic disk position 
angle. This finding, coupled with the result from 
M08 of a similar relationship with respect to the 
disk and nuclear ring ellipticities (fourth panel), 
corroborate s our idea that nuclear rings are cir- 
cular and oriented as their host disk. 

We note three cases that are significantly above 
the scatter and thus do not meet our criteria 
of a round in-plane nuclear ring: NGC 1300, 

NGC 1343, and NGC 5953. The ring in NGC 1300 
is not fully sampled azimuthally (see Figure 2), 
which adds a high level of uncertainty for the pho- 
tometric position angle and ellipticity. NGC 5953 
is nearly face on and morphologically patchy and 
round along the line of sight, a combination that 
creates high position angle uncertainty. 

The residual maps (Fig. 1 and Fig. 2, panel 2) 
also indicate that within the uncertainties of our 
observations (see Section 4.1), the nuclear rings 
are circular. There may be noncircularites at 
a lower level which may be related to the local 
underlying potential (as noted for NGC 4321 - 
Knapen et al. (1995); NGC 1068 - Schinnerer et 
al. (2000); and NGC 5383 - Duval (1977); Sheth 
et al. (2000)). 

5.2. Mass Concentration and Nuclear 
Ring Size 

An indication of the central mass concentra- 
tion can be obtained from the compactness, e.g., 
de Blok et al. (1996). We define compactness as 
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Fig. 5. — Comparisons for ring kinematic vs. ring 
photometric (top), ring photometric vs. disk pho- 
tometric (middle), and ring kinematic vs. disk 
photometric (bottom) position angles. Ellipticity 
comparison of the ring and disk is also plotted. 
The disk position angle is not given in the RC3 for 
NGC 4303, NGC 4314 and NGC 7742. The mean 
error on the disk position angle is unknown, but 
the error increases as the ratio of isophotal major 
to minor axes decreases (logi? 2 s)- We assigned a 
5° error for the smallest ratio in the sample and 
proportionally computed the other errors. Errors 
are smaller than the size of the graphic icons of 
each galaxy. 
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Table 3: Nuclear ring and bar properties. 


NGC 

Nuclear 

RelRing 

Ring 

Qg 

Compactness 

Rt 

Disk 


Activity 

Size 

Width 




sma 




(kpc) 


(10 4 (kms -1 ) 2 kpc -1 ) 

(kpc) 

(kpc) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

473 

- 

0.20 

0.62 

0.12 

1.1 

1.4 

7.4 

1300 

— 

0.02 

0.13 

0.54 

16.5 

0.31 

18.2 

1343 

— 

0.11 

0.31 

0.15 

2.3 

0.3 

11.1 

1530 

— 

0.05 

0.62 

0.73 

8.2 

1.4 

24.5 

2903 

H II 

0.08 

0.01 

0.27 

11.7 

0.1 

16.3 

3351 

H II 

0.04 

0.09 

0.22 

36.1 

0.2 

4.9 

4303 

H II 

0.01 

0.20 

0.28 

23.8 

.2 

14.4 

4314 

LINER 2 

0.05 

0.11 

0.43 

8.6 

0.3 

5.9 

4321 

Transition 2 

0.04 

0.38 

0.22 

15.8 

0.2 

21.7 

5248 

H II 

0.03 

0.49 

0.10 

4.6 

0.5 

20.5 

5953 

LINER 2, Sy2 

0.13 

0.90 

0.1 

9.8 

0.35 

7.7 

6951 

Sy 2 

0.03 

0.23 

0.28 

9.9 

0.25 

13.8 

7742 

Transition 2 

0.18 

0.63 

0.06 

5.8 

0.32 

5.5 


Notes: Identification (col. 1); presence and type of nuclear activity (col. 2, from Ho et al. 1997b when available, 
NGC 1530 and NGC 5953 are from NED); relative nuclear ring size (col. 3) defined as ring semi-major axis (see 
Table 1, col 11) divided by disk semi-major axis (half of the host major disk diameter), Do from NED (col 9); ring 
width (col. 4) defined as the difference between the isophotal inner ring boundary and outer ring boundary, divided 
by the disk semi-major axis; non-axisymmetric torque parameter (col. 5), from Comeron et al. 2009; when no Qb 
value exists, we adopt a known Q g value instead: NGC 1530 (Block et al. 2004) and NGC 5953 (H. Salo, private 
communication); rotation curve compactness (col. 6); rotation curve turnover radius (col. 7); We estimate that on 
average the uncertainty in the compactness is «75kms -1 kpc. Uncertainties in the various size measurements involved 
were also added in quadrature to yield typical uncertainties in the relative ring size and width of ±0.01. Uncertainties 
associated with the non-axisymmetric torque parameter, Q g , were adopted from Comeron et al. (2009) to be ±0.1. 
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Fig. 6. — Left panel: Radial location of the nuclear rings versus the rotation curve turnover radius. The 
locations generally coincide within the uncertainties. The spread is most likely due to the higher uncertainty 
associated with the ridge location (radial ring center) of the wider rings NGC 1343, NGC 4321, NGC 5953, 
and NGC 7742. Right panel: Ring width versus relative ring size. As the relative ring size increases so does 
the ring width. 
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Fig. 7. - Rotation curve compactness versus relative nuclear ring size (panel 1), nuclear ring width (panel 2), 
and the non-axisymmetric torque parameter Q g (panel 3). Panel 4 shows the relative ring size as a function 
of Q g . Typical error bars (refer to Table 3 caption) are within the size of the icons. 


V 2 /Rt, where V is the rotational velocity differ- 
ence between the radial location of the turnover 
radius, R t , and the origin (0,0). The turnover ra- 
dius is defined as the point in the rotation curve 
halfway between the initial steeply rising compo- 
nent and the flat (ter) segment. For the few galax- 
ies with very steeply rising rotation curves, such 
as NGC 2903 and NGC 3351, the turnover ra- 
dius cannot be measured accurately, and hence the 
compactness is uncertain. In these cases, the com- 
pactness should be considered a lower limit. 

Table 3 lists the compactness and turnover lo- 
cations. As expected, the compactness decreases 
with increasing turnover radius (not shown graph- 
ically, but the values are in the Table. In linear 
theory, which is not strictly applicable here but 
yields useful approximations, the turnover radius 
in the rotation curve sets the location of the ILRs, 
and they in turn determine where inflowing gas 
piles up, which ultimately forms the nuclear ring 
(Knapen et al. 1995; Buta & Combes 1996). In- 
deed, the location of the nuclear rings correlates 
well with the turnover radius (see Figure 3 and the 
left panel of Figure 6), which we can identify with 
the phenomenon identified above. 

In panel 1 of Fig. 7 we compare the relative 
ring size, that is the nuclear ring semi-major axis 
divided by the semi-major axis of the host galaxy, 
to the compactness parameter as derived from the 


rotation curve. We find that the relative ring size 
increases as the compactness decreases. 

The relative ring size is compared to the non- 
axisymmetric galactic torque, Q g , in panel four of 
Fig. 7. This torque parameter, originally defined 
by Combes & Sanders (1981) to characterize bar 
strength, quantifies the overall influence of galactic 
non-axisymmetries. Specifically, Q g is the maxi- 
mum tangential force at a given radius divided by 
the radial force. For all galaxies except two, we 
adopt the Q g values from Comeron et al. (2010), 
as noted in Table 3, since they applied a consistent 
computation method for all of the nuclear rings 
seen in our sample. We confirm (Fig. 7) that the 
highest torque values occur in those galaxies with 
the smaller rings relative to the size of their host 
disk (Knapen 2005; Comeron et al. 2010). As the 
torque weakens below Q g = 0.4 any size ring can 
form, but higher torque values (’’stronger bars”) 
progressively limit the nuclear ring size. 

The compactness and the non-axisymmetric 
torque Q g do not correlate (panel 3 of Fig. 7). 
Also, the turnover and Q g do not correlate (not 
shown graphically) . This may be because the non- 
axisymmetric torque is a parameter driven more 
by the azimuthal than the radial distribution of 
mass, whereas the opposite is the case for the 
turnover radius and compactness. 
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5.3. Nuclear Ring Width Comparisons 

The width of a nuclear ring is an important pa- 
rameter because it very much defines how well a 
ring is visible in an Hq- image, and how it is seen 
in contrast with the background emission from the 
host galaxy. Using the ring parameters we have 
now determined, it is interesting to see whether 
the ring width is driven by the physical properties 
of the host galaxy. So, as a first step, we define the 
ring width as the difference between the inner and 
outer ring boundaries, as defined from the emis- 
sion in Ha, and then plot the width of the nuclear 
ring against the relative ring size (right panel of 
Fig. 6). We see a trend that larger nuclear rings 
are wider, which might not have been expected a 
priori. 

We then compare the compactness to the nu- 
clear ring radial width (panel 2 of Fig. 7). There 
is some evidence for a trend of compactness with 
ring width, where a nuclear ring is wider as the 
compactness decreases. We can thus conclude that 
increased compactness leads to both smaller and 
narrower nuclear rings, which is in accord with the 
expectations from linear theory and the location 
of a pair of ILRs near the turnover radius of a 
rotation curve. This clear link between the kine- 
matics of the galaxy and the morphological pa- 
rameters of the nuclear rings is thus further direct 
evidence for a resonant origin of the nuclear rings, 
deeply linked to the underlying structure of the 
host galaxy. This evidence based on kinematics 
supplements earlier evidence based almost exclu- 
sively on morphology, as reviewed and extended 
by Comeron et al. (2010). 

6. Concluding Remarks 

We combined DensePak integral field unit and 
TAURUS Fabry-Perot observations for a sample of 
13 galaxies that contain nuclear rings. Based on 
this sample, we note the following new findings: 

• Nuclear rings are intrinsically circular and 
unaffected by the local non-axisymmetric en- 
vironment that can exist from active star for- 
mation within the rings; they are influenced 
by global phenomenon that can affect the 
potential, such as a bar. 

• The compactness, derived from the rotation 
curves, decreases as the ring width increases. 


• The compactness also decreases as the rela- 
tive ring size increases 

• As the relative nuclear ring size increases so 
does the width of the ring 

• Strong bars with a Qb of 0.4 or greater can 
only be accompanied by small and thin nu- 
clear rings. For smaller torques, the ring size 
can vary. 

We have also reaffirmed several standing find- 
ings: 

• True to resonance theory, the location of 
the nuclear rings correlates well with the 
turnover radius. 

• The highest torque values occur in those 
galaxies with the smaller rings relative to the 
size of their host disk. 

In those cases with sufficient radial coverage, 
we agree that 

• For NGC 4314 residual velocities sharply de- 
crease and plateau within the nuclear rings. 
Ordered velocity patterns increase near the 
outer edge of the ring (both approaching and 
receding) close to both sides of the bar minor 
axis. 

• The residual velocities for NGC 5248 near 
the western side of the bar minor axis reach 
40kms _1 , with a weaker-defined, although 
evident, transition zone as radii approach 
the outer ring boundary. 

• In the case of NGC 1530 large velocity ex- 
cesses exist to the west of the nuclear ring, 
near its exterior edge and then sharply de- 
crease over a small transition zone near the 
outer boundary of the nuclear ring. 

The Fabry-Perot observations were made with 
the William Herschel Telescope operated on the 
island of La Palma by the Newton Group of Tele- 
scopes in the Spanish Observatorio del Roque de 
los Muchaclios of the Institute de Astroffsica de 
Canarias. 

The Integral Field Unit observations were made 
with the WIYN Observatory which is owned and 
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operated by the WIYN Consortium, consisting of 
the University of Wisconsin, Indiana University, 
Yale University, and the National Optical Astron- 
omy Observatory (NOAO). The WIYN is part of 
the Kitt Peak National Observatory located in 
Tucson Arizona. 
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